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ABSTRACT 



Aims. We present predicted Sunyaev-Zeldovich (SZ) properties of known X-ray clusters of galaxies for which gas temperature 
measurements are available. The reference sample was compiled from the BAX database for X-ray clusters. 

Methods. The Sunyaev-Zeldovich signal is predicted according to two different scaling laws for the mass-temperature relation in 
clusters: a standard relation and an evolving relation that reproduces well the evolution of the X-ray temperature distribution function 
in a concordance cosmology. Using a Markov Chain Mote Carlo (MCMC) analysis we examine the values of the recovered parameters 
and their uncertainties. 

Results. The evolving case can be clearly distinguished from the non-evolving case, showing that SZ measurements will indeed be 
efficient in constraining the themal history of the intra-cluster gas. However, significant bias appears in the measured values of the 
evolution parameter for high SZ threshold owing to selection effects. 
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1. Introduction 

Galaxy clusters are the largest gravitationally bound objects 
in the Universe and can therefore provide important infor- 
mation on both the evolution of large-scale structures and 
the astrophysics of collapsed gas and dark matter. In par- 
ticular, the abundance of massive clusters has been shown 
to have the potential to provide strong constraints on the 
cosmological parameters, because it depends exponentially 
on the cosmological growth factor (lOukbir & Blanchard 19921 
IBIanchard & Bartlett 1998l l. Several authors have used clus- 
ter abundance measured through X-ray observations to de- 
rive cosmological constraints, but the results do not always 
agree CBlanchard et al. 20001 [Henry 2004[ IViMiIinin et al. 20091 
IDelsart et al. 201 Oi l. This discrepancy arises because it is difficult 
to find a reliable estimator for the total mass of a cluster, since 
this quantity is not directly observable. Typically, the mass of 
a cluster is estimated with an assumption on the global physics, 
such as a hydrostatic equilibrium of the hot gas in the dark matter 
potential well. This provides a way to connect the total mass at a 
scalable radius to a given observable such as the X-ray temper- 
ature. These indirect measurements together with the results of 



numerical simulations (Ras ia et al. 20 06: Kay etal. 20071 indi- 



cate that at least for local clusters the relation between mass and 
X-ray temperature agree quite well (despite a non-negligible dis- 
persion) with what one would expect from simple scaling argu- 
ments based on the gravitational heating of the gas, with a possi- 
ble uncertainty in its calibration. However, |VaucIair et al. (2003)| 
modeled the abundance of X-ray galaxy clusters at high red- 
shifts (z ~ 0.5) and showed that for observed cluster counts to 
be compatible with a standard cosmology dominated by dark 
energy one would have to consider a significant amount of non- 



standard evolution with redshift in the M - T relation. A simi- 
lar result was obtained recently by Delsart et al. (2010). Indeed, 
it has become generally accepted that non-gravitational heating 
of the gas is significant and can explain the observed departure 
from the L-T scaling relation at z=0 (I Voit 2005T l. This complex 
source of heating could also induce an evolution in the M - T re- 
lation ( Nath & Majumdar 201 1) , which may remain undetected 
even if distant clusters are already in hydrostatic equilibrium. 
Within this reasoning, one can consider that a given scaling law 
can still hold true at higher redshifts if a non-standard evolu- 
tion factor is introduced. This factor could be constrained with 
good observational data up to high redshifts. Of particular inter- 
est is the cross-analysis between measurements of the Sunyaev- 
Zeldovich effect and X-ray temperature, because both data probe 
the physics of the hot gas and are related to the overall cluster 
mass through different scaling relations. 



The 



thermal 



Sunyaev-Zeldovich 



effect 
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(Sunyaev & Zeldovich 1972 1 has become a very promising 
way to observe and detect new galaxy clusters. This effect 
corresponds to a shift observed in the Cosmic Microwave 
Background (CMB) frequency spectrum, which is caused by 
the inverse Compton scattering of the cold CMB photons 
by the high-energy electrons that compose the intramedium 
gas of galaxy clusters. The associated flux decrement can 
be analytically estimated and characterized by the central 
parameter yo or the total integrated flux Y. This provides a 
different way to investigate cluster abundances for cosmological 
constraints (IBarbosa et al. 199"6l Aghanim et al. 1997 1. The 
SZ signal was already measured for some clusters and used 
to study scaling relations between SZ and X-ray observed 
clusters (Morandi et al. 2007). However, to better understand 
and expand these results, it is important to perform a full study 
of the practical implications of assuming different flux limits 
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for the samples as well as having larger samples and smaller 
observational errors when deriving cosmological results. 

The purpose of the present paper is then to investigate the 
capability of present and future experiments in constraining a 
possible evolution with redshift of the Y - T and M - T scaling 
relations, identifying at the same time the eventual presence of 
bias owing to selection effects or other causes. For this we sim- 
ulate the observed SZ integrated fluxes for clusters with known 
X-ray temperatures in two distinct fiducial models with different 
selection criteria. 

The paper is organized as follows. In Sect. 2 we briefly 
present the theory behind the scaling laws for galaxy clusters 
and a model for their evolution. In Sect. 3 we present the cluster 
sample and its simulated SZ observations for different observa- 
tional limits. The results on cosmological constraints, normal- 
ization and redshift evolution are presented in Sect. 4. Finally, 
we present our conclusions in Sect. 5. 



2. Scaling relations for M -T,Y - M and Y-T 

To test the effect of combined measurements on the SZ inte- 
grated flux Y and the X-ray emission-weighted temperature we 
start by considering a general form for the standard scaling law 
between the cluster's mass and temperature. The departure from 
the expected evolution for pure gravitational heating is taken into 
account by introducing an additional term of(l+z)",a being an 
unknown parameter One can write 

T - ArMM'^'/ (^f' /^'''(l + ^)''" keV , (1) 

Am being the enclosed overdensity inside a cluster compared to 
the mean background density of the universe, or equivalently, for 
the same physical radius: 

T = ArMM%' (^^f' lr^\l + zT keV , (2) 

Ac being the enclosed overdensity compared to the critical den- 
sity. We emphasize that both expressions are physically equiva- 
lent and will proceed using Eq. 1 hereafter One can interpret 
these expressions through a physical dependence of the tem- 
perature on the contrast overdensity at a given radius, the total 
mass at this radius (expressed in units of 10 '^M©), and a nor- 
malization constant (Atm), which depends on A. At the virial 
radius, the overdensity is usually computed from the spherical 
top-hat model for gravitational collapse, and for an Einstein- 
de Sitter universe, one has Am,v = - IStt^. An approx- 
imation to the general form of this quantity was provided by 
[Bryan & Norman (1998) | which allows us to determine the virial 
overdensity for any flat ACDM model. In the above expressions, 
the dependence of the background or critical density with red- 
shift is encoded in the (1 +z) and E{z) factors, respectively, with 
E\z) = Q,„(l + zf + Qr(1 + zf + ilA. 

The above M - T scaling relation (Eq. 1) can be used to 
derive a consistent expected relation between the integrated SZ 
flux Y and T. The thermal SZ effect consists of a deformation 
of the CMB spectrum, which can be quantified by the flux (or 
Compton) y parameter along the line of sight: 

y = 7 UeTedl . (3) 



Integrating y over the solid angle Q subtended by the cluster 
introduces the total Y parameter: 

y= r ydo.^'^s^i- r r nj^Adu (4) 

Jil 'neC dg Jo Js 

where da is the cosmological angular distance to the cluster and 
S is the area of the cluster observed in the sky. Y is an observable 
in SZ surveys and represents a unique way to directly measure 
the total thermal energy of the gas content in galaxy clusters. 
One can replace the integrated temperature by the mass-averaged 
gas temperature and write the remaining integral in terms of the 
total mass of the cluster This allows to derive a general scaling 
law from the above expression: 

Y oc TgMgdf = TJgM,„,d-/ , (5) 

Tg being the mass-averaged temperature of the gas. The gas frac- 
tion at the virial radius of clusters is usually expressed as a frac- 
tion of the universal baryonic fraction { fg - T^). Direct mea- 
surements of the gas fraction up to the virial radius are difficult to 
achieve, but numerical simulations seem to agree with a value of 
T between 0.7 and 0.8 (De Boni et al. 20101) . Here we assumed 
a fixed value of 0.8, consistent with a 15% contribution of stars 
to the baryonic component (Roussel et al. 200Tl l. If the gas were 
completely isothermal in clusters, one would have Tg = Tx, with 
Tx being the measured X-ray temperature and Eq.4 could then 
be used to probe the gas mass directly. However, it has been 
shown that clusters are not isothermal, with the temperature de- 
clining in the outer parts (Vikhlinin et al. 2005 ). In this case, 
one can still assume that the gas temperature follows the scal- 
ing law but with a different normalization so that Tg-^Tx, the 
scaling assumption being that ^ is a constant. This normalization 
can be estimated using the observed profiles for gas mass frac- 
tions jSadat et al. 2005l l and temperature (I Vikhlinin et al. 20051 1. 
together with the theoretical NEW profile for the dark matter 
halo. The calculation yields ^ ~ 0.6, a number that is subjected 
to a substantial uncertainty but that could be further constrained 
by observations. 

Using the mass-temperature relation from Eq.l, the above 
expression can be used to derive the full Y -Tx scaling relation, 
which connects two directly observable quantities. Assuming 
that the SZ contribution from regions outside the virial radius 
can be neglected, we have: 

Y = l.816-10-^^A,U(^j ' r] (1 +zri(i-)/,--^-2 .(6) 

If scaling holds, the existence of some evolution in the mass- 
temperature relation and the normalization between the gas and 
the X-ray emission-weighted temperature could be constrained 
using measurements of both Y and Tx for the same clusters. 
Below we address the constraining power of these measurements 
for known X-ray clusters by considering two fiducial cases: the 
standard scaling relation (a - 0) and a significant deviation 
caused by evolution; like |Vauclair et al. (200 3) and Delsart et al. 
(2010) we adopted a = -1 as a fiducial value. As for the temper- 
ature normalization parameter ^, we consider a constant fiducial 
value of 0.57. 

Unless otherwise noted, we performed our study using the 
best cosmological model derived through the COSMOMCQ 
package, using a combination of the WMAP 7-year data 
(IJarosiket al. 20I0i ). SDSS LRG 7 dReid et al. 2010! ) and the 

' http://cosmologist.info/cosmomc/ 
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Fig. 1. Distribution of temperature and redshift for the cluster sample used for this study and the clusters obtained after selection 
cuts for one Monte Carlo simulation. Blue plus signs represent clusters from the original sample that were excluded after the angular 
size cut; green crosses correspond to undetected clusters assuming a detection threshold of LIO"'' arcmin^ (Jeff) or LIO"^ arcmin^ 
(right). Red dots represent detected clusters under the same detection thresholds. This figure remains essentially unchanged in all 
simulations performed in this work. 



Constitution SNIa sample jKessler et al. 20091) . The best-fit 
model obtained is Qm = 0.317, Qa = 0.682, Qfc = 0.0499, 
h = 0.671, where Hq - 100 h km s"' Mpc"'. We assumed a 
fixed value of 8.21 keV for the normalization factor Ajm, cal- 
ibrated through the analysis of the temperature function of ob- 
served clusters (IDelsart et al. 2010t . 



3. Simulated SZ observations of real clusters in two 
evolution scenarios 

To build the largest possible sample, we made use of the X-ray 
Galaxy Cluster Database (BAX). We selected a sample contain- 
ing all the identified clusters and groups for which temperatures 
are published in the literature with a relative measurement error 
of less than 15%. This makes the sample far from being uniform 
but it serves the purpose of obtaining the largest sample of po- 
tentially observable clusters in upcoming SZ surveys with mea- 
sured X-ray temperatures. Under these conditions, we were able 
to select a total of 438 clusters, with redshifts between 0.003 and 
1. In FiglUwe plot the distribution of temperature versus redshift 
for the global sample, allowing for a better visualization of the 
sample measured quantities. In the same figure we show clusters 
that passed selection cuts described below in the paper 

To perform our study, we started by generating "real" tem- 
peratures for each cluster through a simple Gaussian sampling 
based on the measured mean values and uncertainties. These real 
temperatures were then used as a starting point to derive other 
observable quantities of the sample, based on the scaling laws 
presented above. 

One important feature when measuring the SZ flux is the an- 
gular size of the cluster in the sky. Large clusters will be fully 
resolved and an exact measurement of the Y parameter up to the 
virial radius will require a good modeling of the cluster internal 
structure and pressure profile. On the other hand, smaller clusters 
will become nearly or completely unresolved and the virial mea- 
sured integrated flux will thus become less model-dependent. 

In Fig.|2]we show the distribution of the angular radius for 
the clusters in our sample (for one realization), which presents 



a peak at about 8 arcmin and extends to values as high as 50 
arcmin. These values were derived from Eq. [T]together with the 
condition My - 4/3npi„6vRy and the conversion to angular size 
using the cosmological angular distance. One can see that by 
introducing a significant evolution in the M - T relation, the 
overall size distribution is still preserved. We decided to only 
select clusters with an apparent radial size of less than 10 arcmin, 
a value about twice higher than the expected angular beam size 
for the Planck SZ survey Q. This choice was made because once 
we used the Planck resolution as a sample limit, the number of 
clusters would be drastically reduced to 30. Adopting this nearly 
resolved margin, and using 15 different realizations, the sample 
after this cut in angular size was composed of 199+1 clusters. 
This low variance reflects the good precision in the measured 
temperatures of the initial cluster sample. 

We proceeded then to construct a "mock" catalog of the ob- 
served Y values. Because we are interested in the dependence 
of the Y - T relation on redshift, which is modeled as a power 
law, we chose to work in logarithmic units for both Y and T. 
In this way, the "real" \og{Y) was obtained by taking the loga- 
rithm of Eq. [T] and introducing the cosmological parameters in 
our fiducial model together with the "real" temperature of each 
clusters. Furthermore, we modeled an intrinsic dispersion of the 
scaling law of 10% by Gaussian sampling around this first esti- 
mate of log(y). This value corresponds to the typical result for 
the dispersion in M - Tx or M - F obtained in numerical simula- 
tions jKravtsov et al. 20 06^. Moreover, the recent analysis from 
the Planck collaboration (Planck Collaboration 201 lai for low- 
redhift clusters seems to support a value of 12% at the Icr level, 
consistent with 10% at the 2cr level. 

To the values obtained after the introduction of intrinsic 
scatter we added random measurement errors based on the 
estimated flux detection threshold for the Planck SZ survey. 
Following |Mehn et al. (2006)| and the recent Planck data release 
(IPlanck Collaboration 20 1 IcV we considered the value for the 
threshold to be of 1.10"^ arcmin^. Moreover, we assumed the 



http://www.rssd.esa.int/index.php?project=PLANCK&page=perf_top 
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Fig. 2. Histogram of the angular size of our cluster sample as 
expected from standard M -T relations (top) and accounting for 
a strong evolution in M - T with a = -1 {bottom). 
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Fig. 3. Histogram of the Y values in logarithmic form obtained 
in one simulation, using the scaling law defined in Eq. |6] with 
a = (top) and a = -1 (bottom). An intrinsic dispersion of 10% 
is assumed. 



measurement errors to have a mean value corresponding to 20% 
of the detection threshold, which corresponds to a signal-to- 
noise ratio of 5. After this procedure we obtained a simulation 
of the observed integrated flux parameters using our two fiducial 
models for Y -T and M -T, with a distribution (for one realiza- 
tion) presented in Fig. 3. This distribution presents a peak around 
10"^ arcmin^ and extends up to 4.10"^ arcmin^, thus confirming 
the importance of having a low detection threshold. In fact, by 
setting the detection threshold at 1.10"^ arcmin-, the total num- 
ber of detectable clusters is reduced to 79+4 (estimated from 15 
realizations). If o- = -1, the sample size increases to 107+4. 



4. Constraining power of F - r observations 

In this section we explore the use of data on X-ray temperature 
and SZ flux to derive constraints on both the cosmological pa- 
rameters and the parameters of the scaling laws in clusters. The 
"mock" Y data derived in the previous section were used to iden- 
tify the constraining power and possible presence of bias when 
combining these two information sources. 

We are thus interested in fitting the data with a scaling law 
(Eq|6]l, which in logarithmic form can be written as 



iogr = ^: + iog^+^iogr-^(i+a)iog(i+z) 



1 



■ log Q„ - log h-2 log t/a(Q,„, Qa, /!) . 



(7) 



where is a constant factor accounting for terms like fgas, phys- 
ical constants and the normalization Ajm, while the variation in 
these parameters as well as the SZ contribution beyond the virial 
radius can be absorbed in our parameter ^. The other terms are 
explicitly shown to better illustrate the linear dependence on the 
observables log T, log(l + z). We considered two data sets, one 
constructed from a standard scaUng law (a - 0), and a second 
data set comprising a strong non-standard evolution {a - -1). 

An important aspect to take into account when fitting this 
model is that both dependent and independent data present non 
negligible measurement errors. In this case, Y, T, and z are mea- 
sured with some degree of uncertainty and, although the cluster 
redshifts z are usually measured with high precision, this is not 
the case for the temperature and SZ flux. It is necessary then 
to use a statistical estimator that takes this into account by con- 
structing an effective using a total least squares regression 
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Table 1. Table summarizing the mean posterior and 68% confidence intervals of the constrained parameters. 
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t 0.05(±0.04) 


0.54 ± 0.04(±0.03) 


0.56 ± 0.04(+0.04) 


0.57 ± 0.02(±0.03) 


a 


-0.27 


±0.21(±0.18) 


-1.29±0.15(±0.14) 


-0.04±0.09(±0.10) 


-0.99 + 0.09(±0.11) 


Fiducial model 




a = 


or = -1 
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a = -1 


Threshold (arcmin^) 




1.10^3 


1.10^3 


1.10"^ 


I.IO^'' 


O" disp 


0.078 ± 0.040 


0.068 ± 0.049 


0.113 ±0.018 


0.103 ±0.019 




0.56 H 


t 0.03(±0.05) 


0.55 ± 0.04(±0.04) 


0.57 ± 0.04(+0.03) 


0.57 ± 0.04(±0.03) 


a 


-0.22 


+ 0.23(±0.18) 


-1.29±0.17(±0.17) 


-0.02 + 0.09(±0.10) 


-0.98 ± 0.09(±0.09) 



Notes. Table summarizing the obtained constraints on scaling parameters and intrinsic dispersion. Four combinations of fiducial parameters and 
SZ detection limits are presented, using two approaches to the intrinsic dispersion of the scaling laws. The upper part of the table shows the results 
considering the same fiducial value as in the simulations (cTrf/jp = 0.1) . The lower part shows the results when cr^;,,, is derived from the simulated 
sample itself. For the parameters a and f the first quoted errors correspond to the mean value of the 68% confidence errors derived from each 
realization; the second quoted uncertainties represent the standard deviation of the mean central values from the fifteen realizations. The quoted 
uncertainties given for crjisp are the standard deviation derived from the fifteen realizations. 



technique. If F,, T,, and zi are the measured values correspond- 
ing to the same cluster (F, corresponding to our simulated val- 
ues) and F(Ti, Zi, ^, a) is the linear model as on the right hand 
side of Eq|7] the total least-squares can be written as 



i crt, -H 



■F(Ti,Zi,^,a)Y 



{(J 



1,1 disp 



(8) 



This expression assumes that the measurement errors cry, and 
(TT,i are uncorrected and neglects the contribution from redshift 
measurement eiTors (cr- ,). Moreover, it takes into account the 
intrinsic scatter (crj,-,^,) around the scaling law. 

The above estimator was used to constrain a flat ACDM 
model with six free parameters (ilcoMh^, ^th^, h, t, n,. A,) 
plus the two free parameters of our scaling model ^ and a. 
We added a module to the public Markov Chain Monte Carlo 
code COSMOMC, and performed a joint constrain using CMB, 
LSS, SNIa and our SZ-Xray data. This combination can pro- 
vide strong constraints on cosmological parameters required to 
reduce parameter degeneracies, which improves the constraints 
on cluster-related parameters. 

As stated above, the nature of Monte-Carlo simulations is 
prone to an intrinsic variability when drawing different realiza- 
tions. For this reason, we computed a total of 15 realizations for 
each considered fiducial model and detection threshold and used 
them to perform the MCMC chains. The mean values derived 
from the posterior distribution on the parameters of intere st were 
then averaged to account for the intrinsic variability. For consis- 
tency, we compared the standard deviation for each parameter 
over the 15 realizations and verified that it was consistent with 
the typical Icr eiTor for a single realization. 

Another important aspect when using this kind of data is the 
estimation of the intrinsic dispersion o-^isp to be accounted for 
in Eq. [8] Similarly to what is done with the SNIa Hubble di- 
agram (Kowalski et al. 2008; Hicken et al. 2009), we estimated 
this quantity using a goodness-of-fit criterion and residual anal- 
ysis. Although we knew by construction the real value for crj,.5p, 
we adopted an iterative method to find the intrinsic dispersion 
using a maximum-likelihood criterion. In this case, the "correct" 
value should verify xlff ~ N - p, with being the total number 
of clusters and p the number of fitted parameters. To implement 
this, we ran MCMC chains up to 10000 points for a given data 
set and assuming an initial estimate of cr^/jp = 0. After having a 
reasonable estimate of the minimum value and the associated 
parameters, the residuals were fitted to obtain a first estimate of 



the intrinsic dispersion, which was then used to repeat the pro- 
cedure. After a small number of iterations we let the MCMC 
run until full convergence, which produced a best-fit model with 
X^ X N - p. This allowed us to study the impact of the assumed 
intrinsic dispersion in the constrained parameters. 

Table 1 summarizes the obtained results for the constrained 
parameters in the two fiducial models, assuming different values 
for the flux detection threshold. These values account for the av- 
erages and standard deviations obtained from the 15 realizations 
made for each case. We explicitly show the differences obtained 
when the correct intrinsic dispersion of 0.1 was assumed (up- 
per part of the table) or when it was estimated using the method 
described above (lower part). In the later case, the values esti- 
mated for cr^iisp were also averaged over the 15 realizations and 
we provide the obtained mean and stdev values. 

The results show a noticeable dependence of the scaling pa- 
rameter a on the assumed detection threshold. When the ini- 
tial value of 1.10"^arcmin^ is assumed, the marginalized mean 
values for the evolution parameters are slightly underestimated, 
with the fiducial values (used to construct the data set) lying out- 
side the derived 68% confidence intervals. This is most likely be- 
cause the clusters with lower values of Y (which are not detected 
at this threshold) are often those at higher redshifts or those pre- 
senting a significant deviation from the mean value. When a 
higher threshold value is considered, the absence of these un- 
detected clusters will cause the slope of the Y - z relation to be 
tilted towards lower values. A similar but smaller effect occurs 
with the normalization parameter 

To test this explanation, we used a lower value for the thresh- 
old in both cases to check if we could retrieve better estimates 
of the input fiducial parameters. Indeed, reducing the threshold 
to 1.10"^ arcmin^ causes the number of detected clusters to in- 
crease significantly (detecting all the nearly unresolved clusters 
in our initial sample), as shown in Fig.l. This causes the con- 
straints on the scaling parameters to become essentially unbiased 
and with a precision on a improved by a factor of 2 to 3, depend- 
ing on the evolution scenario. In Fig. 4 we present the 2D confi- 
dence regions for four realizations with different fiducial values 
of a and detection threshold. The improvement in the constraints 
becomes evident when assuming a lower detection threshold. We 
note that the threshold does not have to be reduced to a value as 
low as 1.10"^ arcmin^ to lead to a significant increase in the to- 
tal number of detected clusters. Indeed, in Fig. |3]one can see that 
only a small number of clusters remain undetected if one consid- 
ers the threshold to be 1 . 10 arcmin^. However, the few clusters 
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Fig. 4. Two-dimensional marginalized constraints on the scaling 
model parameter, representing 68% confidence regions for two- 
parameters. Blue (upper) contours were derived from data simu- 
lated with a - 0, and red (lower) contours using a - -\. Filled 
contours assumed a Y detection threshold of 1.10"^ arcmin-, 
while dashed curves assumed this value to be 1.10"^' arcmin-. 



that remain undetected have intermediate to high redshifts (be- 
tween 0.5 and 1), which can be important to constrain evolution 
effects on scaling laws. Also, a lower threshold ensures better 
precision in individual cluster measurements. 

Nevertheless, we computed MCMC chains considering a 
1.10"^ arcmin^ threshold and obtained a fairly intermediate re- 
sult, i.e., larger uncertainties than those obtained with a 1.10"'' 
arcmin^ threshold, but nearly consistent with the assumed fidu- 
cial values at Icr level. 

One can infer from Table 1 that the method used to esti- 
mate the intrinsic scatter of the scaling law (cr^,.5p) is able to 
retrieve the fiducial value within the Icr statistical variability of 
the Monte Carlo simulation for all studied cases. The decreasing 
of the detection threshold and the related inclusion of more clus- 
ters allows us to improve the accuracy with which the fiducial 
value is retrieved. We also note that the obtained constraints are 
highly stable, independent of the assumed method to account for 
the intrinsic dispersion (Jdisp- 

5. Conclusions 

We performed a study of the capacity of using measurements 
on the SZ integrated flux to constrain a possible non-standard 
evolution of the scaling relations of galaxy clusters. Based on 
simulated data derived from measured X-ray temperatures, we 
showed that a sample of 199±1 clusters with measured temper- 
atures (comprising 79+4 detectable and nearly unresolved clus- 
ters within Planck's SZ observational "settings") can indeed pro- 
vide meaningful information on this subject, and distinguish be- 
tween a strong non-standard evolution and the standard gravita- 
tionally driven scaling law. 

We showed that even if some intrinsic dispersion is present 
in the actual scaling law for Y - T, its amplitude can be esti- 
mated and extracted from the data with an acceptable precision, 
using a goodness-of-fit analysis. More important is that this in- 
trinsic dispersion does not affect the parameter estimates, allow- 



ing for a possible distinction between standard scaling models 
and strongly evolving ones. 

Using a conservative value for the SZ flux measurement 
threshold (1.10^^ arcmin^), the relative uncertainty obtained on 
the evolution index a is significant, because of the low number of 
unresolved detectable clusters. Furthermore, we showed that by 
imposing this relatively high flux limit, the obtained constraints 
on this parameter present a significant bias that needs to be taken 
into account. This implies that although the Planck SZ survey 
can provide an important distinction between non-evolving and 
strongly evolving Y - T relation, its results can be affected by a 
non-negligible systematic error Precision temperature measure- 
ments through an X-ray follow-up of newly SZ discovered clus- 
ters could in principle help to reduce parameter uncertainties, 
and we showed that lowering the Y detection flux limit allows 
one to obtain unbiased estimates for the evolution parameters. 

After the submission of this paper, the Planck collabo- 
ration has released its early results in a set of papers, with 
several ones focusing on the SZ detection and measure- 
ments of cluster properties. (Planck Collaboration 201 laj 
Plan ck Collaboration 201 Ibl Planck Collaboration 201 IcV 
Although these comprise an analysis of the scaling relations at 
low redshift, the method used to treat the Planck data is based 
on the assumption of a universal profile and scaling law with 
parameters derived from X-ray data alone. In these analyses, the 
measurements are restricted to the central part of clusters . Our 
work shows that it should be possible to address the possible 
non-standard history of the thermal content of clusters, provided 
that the bias induced in Y by selection effects is properly taken 
into account. 

In the wake of a new era for SZ cluster studies, the use of 
SZ flux measurements together with X-ray data can undoubt- 
edly shed some light on the evolution of the properties of these 
massive objects, which can result in important consequences for 
cosmological applications. 
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